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Abstract 

We argue that with an increase of the collision energy, elastic photoproduction of p mesons on 
nuclei becomes affected by the significant cross section of photon inelastic diffraction into large 
masses, which results in the sizable inelastic nuclear shadowing correction to (j 1 a-^ p a and the re¬ 
duced effective p-nucleon cross section. We take these effects into account by combining the vector 
meson dominance model, which we upgrade to include the contribution of high-mass fluctuations 
of the photon according to QCD constraints, and the Gribov-Glauber approximation for nuclear 
shadowing, where the inelastic nuclear shadowing is included by means of cross section fluctuations. 
The resulting approach allows us to successfully describe the data on elastic p photoproduction on 
nuclei in heavy ion UPCs in the 7 GeV < W JP < 46 GeV energy range and to predict the value of 
the cross section of coherent p photoproduction in Pb-Pb UPCs at y' s^n = 5.02 TeV in Run 2 at 
the LHC, dapbpb^. p pbpb(y = 0 )/dy = 560 ± 25 mb. 
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1. Introduction 

At high photon beam energies E y , the photon participates in the strong interaction with hadrons 
through its fluctuation into hadronic components. The lifetime of these fluctuations is character¬ 
ized by the coherence length l c = 2E 1 /M 1 2 , where M is the mass of a given hadronic component. 
With an increase of E 1 , l c increases and becomes larger than the target size Rp for progressively 
heavier hadronic fluctuations of the photon, which is manifested in an increase of photon inelas¬ 
tic diffraction into large masses. This means that the photon can be represented as a coherent 
superposition of hadronic fluctuations interacting with the target with a wide spectrum of cross 
sections. This picture can be implemented in terms of either hadronic (in particular, vector mesons) 
or quark-gluon degrees of freedom. 

In the 60’s and early 70’s, it was suggested that the observed hadron-like behavior of a photon 
in photon-hadron interactions can be represented by the vector meson dominance model (VMD), 
which assumes that the photon fluctuates into p, uj and 0 mesons that subsequently interact with 
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hadrons 0,3- The VMD model successfully explained the behavior of the pion form factor, 
certain features of the nucleon form factors at small momentum transfers and the major part 
(fa 80%) of the real photon-nucleon total cross section a lp . Combining VMD, the constituent quark 
model [3] and the Regge-Gribov theory of high energy hadron-hadron scattering, the soft Pomeron 
phenomenology had been developed and actively used [ 3 , 0 to study properties of the vector meson 
interaction with nucleons in light vector meson photoproduction and electroproduction at small 
photon virtualities Q 2 . One of the key features of this approach is the assumption that a p n = cr-n-N, 
which is based on the additive quark model (cr pN and oyjv are the total p-nucleon and pion-nucleon 
cross sections, respectively). In a wide range of pion energies, the total pion-nucleon cross section 
is described well by a sum of the soft Pomeron and the secondary Reggeon exchanges 00 (we 
refer to this model as DL94). 

An increase of the photon virtuality Q 2 leads to a gradual transition from the soft nonpertur- 
bative QCD regime to the perturbative one, which is clearly revealed in HERA measurements of 
vector meson electroproduction on the proton at high energies, for the review and references, see, 
e.g. [bj. To explain in a consistent way the fa 20% discrepancy between the experimental value 
of cr 7P and VMD predictions and the behavior of photon-induced processes with an increase of 
Q 2 , new theoretical approaches have been developed. On the one hand, within the framework of 
hadronic description, the VMD model was generalized on the basis of the mass dispersion relation 
to include higher-mass resonances with their diagonal and non-diagonal transitions 0, iiQ; 
the resulting approach is referred to as the generalized vector meson dominance (GVMD) model. 
On the other hand, in the QCD framework, the photon wave function in the strong interaction 
can be modeled as superposition of quark-antiquark pairs (dipoles) and higher Fock states, which 
interact with the target. The approach in which one treats the photon as a quark-antiquark pair 
is called the color dipole model (CDM). 

From the point of view of the quark-hadron duality, these two approaches should be in principle 
equivalent, but this equivalence is destroyed in their practical realization. In particular, to apply 
the CDM in the nonperturbative domain, for example, for the description of photoproduction of 
light vector mesons, the approach should be generalized to take into account more complicated 
than qq states of large masses and also to model the cross section of the interaction of large-size 
dipoles with nucleons. On the other hand, the GVMD model includes coupling constants of the 
photon to higher-mass resonances and amplitudes of their diagonal and non-diagonal transitions, 
which are in general unknown. As a result, both approaches require to use phenomenology and 
engage additional experimental information. 

Coherent photoproduction of light vector mesons on nuclear targets at low and intermediate en¬ 
ergies has been considered as theoretically well understood within the framework of the VMD model 
and the Glauber theory of multiple scattering 0 taking into account the finite coherence length 
(for brevity, we refer to the resulting approach as VMD-GM). Recently the ALICE collaboration 
presented results on exclusive p meson production at the central rapidity in Pb-Pb ultraperipheral 
collisions (UPCs) at y/s^N = 2.76 TeV [l2[. The value of the cross section of coherent p meson 
photoproduction on lead at the photon-nucleon center-of-mass energy of W rp fa 46 GeV extracted 
from this measurement is a 1 pb^-pPb ~ 2 mb, which is very close to the 'yAu —y pAu cross sections in 
the range of energies of W rp fa 7 —12 GeV obtained by the STAR collaboration in ion-ion UPCs at 


RHIC 13, 14, [15]. While the calculations of the standard elastic nuclear shadowing in the Glauber 
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model capture the bulk of the nuclear suppression by the factor of approximately four, the ALICE 
experimental values are still significantly, by the factor of approximately 1.5 — 1.7, smaller than 
the VMD-GM predictions of [3, H3, |l8| . It is also important to emphasize that the VMD model 
with the standard values of the photon-p meson coupling constant f p and the p-nucleon cross 
section cr p ^ overestimates the most recent HI data on diffractive p photoproduction on the proton 


at HERA 19] by the factor of approximately 1.3. 


The aim of the present paper is to modify both VMD and Glauber models at high photon 
energies to take into account an increasing role of hadronic fluctuations of the photon interacting 
with different strength and having a wide range of masses. 

At high energies, the Glauber model is substituted by the Gribov-Glauber approach [20J, which 
takes into account both elastic and inelastic diffraction in the intermediate states contributing to 
the shadowing correction. Note that in spite of very different pictures of the space-time evolution of 
the scattering process at moderate and high energies, the expressions for the shadowing correction 
to hadron-nucleus cross sections look similar in the two approaches. In particular, the nuclear 
shadowing is calculable in terms of the elementary hadron-nucleon diffractive cross section, which 
includes the clastic hadron-nucleon scattering leading to the elastic nuclear shadowing correction 
of the Glauber model and the inelastic hadron-nucleon scattering. With an increase of energy, 
inelastic diffraction becomes essential and this results in Gribov’s inelastic nuclear shadowing 
(GINS) correction. 

While an increase of the GINS correction with energy is well-known in the discussed energy 
range, it is not often appreciated that the relative magnitude of the effect is larger for projectiles 
with smaller hadron-nucleon cross sections. Indeed, in proton-nucleus scattering, the GINS cor¬ 
rection to the total proton-nucleus cross sections is found to be small 21, 22|. As a consequence, 


the Glauber model is now widely used in the analysis and interpretation of the data on high en¬ 
ergy heavy ion collisions at RHIC and the LHC. For pion-nucleus scattering, the analysis of total 


pion-deuteron 23] and pion-nucleus 24] cross sections demonstrated that the inelastic nuclear 


shadowing correction is larger than that in the proton-nucleus interaction. We argue that in the 
high energy p meson photoproduction on nuclear targets the GINS correction can be even stronger 
than in the pion-nucleus case. 

We pointed out above that an increase of the photon inelastic diffraction into large masses 
can be understood as resulting from photon fluctuations with the large invariant masses. These 
fluctuations lead to a certain decrease of the yN —$■ pN cross section compared to the expectations 
based on the VMD and the additive quark model assumption cr pN = g^n- Indeed, in the GVMD 
model, the partial cancellation between the diagonal pN —y pN and the non-diagonal p'N —y pN 
transitions, which is required to restore the Bjorken scaling of the total virtual photon-proton 
cross section, naturally leads to a decrease of a lN ^ pN . In the CDM, the quark-antiquark dipoles 
with the large relative transverse momentum k t and the longitudinal momentum fraction z ^ 0,1 
are characterized by the large invariant mass, the small transverse size d t and, hence, the reduced 
cross section due to color transparency of QCD. Since the overlap integral between the photon 
and p meson wave functions has a stronger support at small dipole sizes dt than the square the 
p meson (pion) wave function due to the point-like photon coupling to quarks, the contribution 
of such dipoles somewhat decreases In the framework of the VMD model, this reduction 

can be explained either by an increase of f p or a decrease of g p jy. Since the former is very well 
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constrained by the measured p —> e + e~ decay width, one concludes that <j p n should be somewhat 
reduced, cj p n < ct^n- 

This results in an overall decrease of the ylV —> pN cross section and the relative increase of 
the GINS correction. Besides, the effect of the inelastic shadowing correction in the elastic cross 
section is larger than that in the total cross section. Hence, these modifications of the VMD-GM 
approach noticeably reduce the cross section of coherent p photoproduction on nuclei. 

To implement these effects in our calculations, we (i) use the framework of cross section fluctu¬ 


ations [25|, [26, |27|, 128[, |29|, [30] taking into account quark counting rules for the probability of cross 


section fluctuations for small a and the information on inelastic diffraction in photon scattering off 
a nucleon target, and (ii) modify the VMD model to effectively include the effect of the reduction 
of the p-nucleon cross section. The resulting approach allows us to describe well the data on elastic 
p photoproduction on nuclei in heavy ion UPCs in the 7 GeV < W lp < 46 GeV energy range in a 
way consistent with the yp —>- pp HERA 2006 data [lj|, and, thus, to explain away the discrepancy 
between the data on the yA —> pA cross section and its theoretical description in the VMD-GM 
approach. 


2. Cross section of p photoproduction off nuclear targets from STAR and ALICE 
UPC measurements 


The basic expression for the cross section of vector meson V photoproduction in nucleus-nucleus 
UPCs (for review of the physics of ultraperipheral collisions and references, see 31]) reads: 


(>n ' ' ^ ^ = N l/A (y)a^ A ^vA(y) + N l/A (-y)<J lA ^ VA {~y ), ( 1 ) 

where y = ln(2 uj/My) is the rapidity of the vector meson; oj is the photon energy; My is the vector 
meson mass; a lA ^y A is the cross section of exclusive coherent photoproduction of V on nucleus 
A; N 1 / A (y ) is the photon flux. The presence of two terms in Eq. (Jl]) is the reflection of the fact 
that each colliding nucleus can serve as a photon source and as a target. 

The flux of photons produced by a fast-moving point-like charge is well-known from classical 
electrodynamics: 


AT ( \ 2 ^ 2 «e.m. 

^ 7 /AW = 


7T 


Cfto(C)A-i(C) - j {KUO - Kq(O) 


( 2 ) 


where Z is the nucleus charge; a e . m . is the fine-structure constant; K 0 and K\ are modified Bessel 
functions of the second kind; ( = Myb m \ n e~ v /(2 07 ,), where y l is the nucleus Lorentz factor and 
b m hi is the minimal impact parameter of the nucleus-nucleus ultraperipheral collision. The photon 
flux calculated using Eq. (J2]) with 6 min = 2R A ( R A is the nuclear radius) reproduces with the 
precision of a few percent a more accurate calculation, which takes into account the nuclear form 
factor and the suppression of the strong nucleus-nucleus interaction calculated using the Glauber 


model, see the discussion and references in [32j, |33]. In the latter work, the very good precision of 


the theoretical calculation of N^/A^y) was explicitly demonstrated by the measurement of neutron 
emission in electromagnetic dissociation of Pb nuclei at the LHC at y's/viv = 2.76 TeV. 
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At y — 0, the two terms in Eq. (HD are equal and one can then unambiguously determine the 
7 A —> pA cross section using the experimental values of the da AA ^. AAp (y = 0 )/dy cross section 
measured by the STAR collaboration at RHIC 113, 14, 151 and the ALICE collaboration at the 
LHC 0: 

1 da AA^AApiy = 0 ) 


^*7 a —> p /i (Hiyp) 


2 N l/A (y = 0 ) 


dy 


( 3 ) 


where W 7P = W 7P (y = 0) = \/2M p mNyL with M p = 770 MeV being the p meson mass and uln 
the nucleon mass. The qA —> pA cross section determined this way is presented in Table H] as a 
function of the corresponding W 1P . Note that similarly to the invariant energy of AA collisions 
per nucleon ,/snn, throughout our paper W lp denotes the invariant photon-nucleus energy per 
nucleon. 


Nuclear target 

W 7P: GeV 

^7 A-^fpA'i m b 

qAw —> pAu 

6.96 

2.08 ±0.33 

q Au -A pAu 

10.04 

1.9 ±0.76 

q Au -A pAu 

12.46 

1.58 ±0.25 

q Pb -A- pPb 

46.28 

1.97 ±0.21 


Table 1: The cross sections of p pliotoproduction on nuclear targets extracted from the STAR 
ALICE UPC measurements [12J. 


SHE! and the 


3. Nuclear shadowing in p photoproduction on nuclear targets 

In the Glauber model, the cross section of coherent p photoproduction on nuclei reads jH]: 


a. 


VMD 
7v4 —ypA 


da. 


ip->pp 


(t = 0) 


dt 


dt 


d 2 be^ [ dzp A (b, z)e^ z e- {1 - ill)Z ^ ^ dz ' pA ^ z,) 


( 4 ) 


where a p n is the total p-nucleon cross section; q is the ratio of the real to the imaginary parts of 
the p-nucleon scattering amplitude; q± and q\\ are the transverse and longitudinal components of 
the momentum transfer to the nucleus, respectively; p A is the nuclear density normalized by the 
relation f dz d 2 bp A (b, z) = A. The nuclear density is well known from studies of elastic electron 


and proton scattering on nuclei at intermediate energies (see, for example, [34|, [35|, (36]). In our 


calculations we used the Hartree-Fock-Skyrme nuclear density which describes the root-mean- 
square radii of nuclei across the periodic table with the precision better than 2% 37]. The minimal 


where W rp is the invariant 


momentum transferred squared is f min = —(M 2 7tin /W 2 p ) 2 = —( 911 ) 
photon-nucleon energy (W 2 p = 2 m,]yE 7 + m 2 N in the laboratory frame). In the VMD model, the 
forward elementary qp —>■ pp cross section in Eq. (J1D can be related to the total pN cross section, 
a pN , using the optical theorem: 


da. 


7 p^rpp 


(* = 0 ) 


dt 


16n 


j) (l + h>>, 


( 5 ) 
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where f p is the 7 — p coupling constant (f 2 /Air = 2.01 ± 0.1) fixed by the T(p —> e + e ~) width of 
the p —> e + e~ decay. Note also that the effect of t m ; n 7 ^ 0 in da lp ^ pp (t = 0 )/dt is negligibly small 
and can be neglected. 

At high photon energies, neglecting the effects of q\\ 7 ^ 0 and rj_ 7 ^ 0 in Eq. ((4]) (in the 5 GeV < 
W 1P < 50 GeV range, \rj\ < 0.1 and thus can be safely neglected [4j]), one can write it in the form 
of the optical limit of the Glauber model: 


VMD _ 
a 'yA^pA ~ 



d 2 b 


o -Z£*T A (b) 



( 6 ) 


where Ta(5) = f dzpAip , z). 

The interpretation of Eq. (J6]) is straightforward and well-known: the incoming photon trans¬ 
forms into a p meson far before the target and then the p meson coherently interacts with the 
nucleons along its trajectory. The Glauber model takes into account the elastic nuclear shadowing 
effect, which depends on the total pN cross section. Based on the additive quark model, it is 
generally assumed that ct p n(W) = ovv(IU) = [c j n + N (W ) + cr 7r -j V (kE)]/2, where W is the invariant 
p meson (pion)-nucleon energy. The experimental pion-nucleon cross sections in a wide range 
of energies are well described as a simple sum of the soft Pomeron and the secondary Reggeon 
exchanges [ 4 J. Correspondingly, this resulted in a simple form (DL94) for the energy dependence 
of a p n |5(]: 

a pN (W) = 13.6 w 2(ap(0)-1) + 31.8 W 2( “ i?(0) - 1) , (7) 


where or(£) = 0.55 + 0.93t is the Reggeon exchange trajectory; ap(t ) = 1.0808 + 0.25f is the 
soft Pomeron trajectory characterizing the high energy behavior of soft hadron-nucleon processes. 
Since the DL94 model did not fit well the forward p photoproduction cross section on the proton 
measured to that time in low energy experiments [38|, |39|, 


40] and in the ZEUS experiment 4ll |42 
43] at high energies, Donnachie and Landshoff suggested [44[ to simply renormalize the forward 
cross section [Eq. (J5])] by the factor of 0.84 motivated by possible corrections to the 7 — p coupling 


constant (this renormalization was consistent with large experimental uncertainties of the data). 

In Fig. [1] we compare the total cross sections of coherent p photoproduction extracted from the 
STAR and ALICE measurements (see Table [T]) to those calculated in the impulse approximation 
(IA), when all nuclear effects except for coherence are neglected (blue dot-dashed line), and in the 
Glauber model (red dashed line) using the DL94 pN total cross section (Eq. (TT|) ). From the figure 
one finds that the VMD-GM with the DL94 model for ct p n predicts the suppression of o^A-^pA 
by approximately a factor of four compared to the IA calculation, but it still overestimates the 
experimental cross sections by the factor of 1.5 — 2. Besides, the energy dependence is different: 
while the calculated cross sections slowly grow with energy, the experimental values slightly de¬ 
crease or stay almost constant. Note that the calculated values of the 7 Au -A- pAu cross section 
are smaller than those for the lead target by approximately 5% for all energies. Hence, we neglect 
this difference throughout our paper and perform our calculations for lead keeping in mind the 5% 
reduction of the nuclear cross section when we compare our calculations with the STAR data. 

To check the accuracy of the Glauber model calculations in Eq. ([HD in combination with the 
DL94 pion-nucleon cross section, we calculated the hadron-nucleus total and inelastic cross sec- 
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Figure 1: The yA —» pA cross section as a function of W 7P . The VMD-GM (red dashed curve) and VMD-IA 
(blue dot-dashed line) predictions for a 208 Pb target based on the DL94 parametrization of the pN cross section 
are compared to the experimental values extracted from the STAR and ALICE UPC measurements. 


- 7+Pb—>p+Pb in VMD-GM-DL94 

- - 7+Pb—>p+Pb in VMD-IA-DL94 



• STAR 7+Au^p+Au 1 
a ALICE 7+Pb^p+Pb ) 


20 40 60 80 

W™ GeV 


tions for the neutron and pion projectiles in the Glauber approach: 


<1 = 2/^6 


1 _ e -^FT A {b) 


a hA ~ 


d 2 b[ 1 


^—^knTa{V) 


]• 


( 8 ) 


The neutron-nucleon cross section a n N is estimated using the additive quark model counting rule 
relation [3j a n n = 3/2cr nN , where the pion-nucleon cross section is given by Eq. (J7J). The results 
of our calculations are compared to the data 24, |45|, |4b|, [47] in Fig. [21 One can see from the figure 
that the calculations agree very well with the measurements. This means that the reasons of the 
disagreement of similar calculations of the 7 A —> pA cross section with the STAR and ALICE 
data are in specifics of the light vector meson photoproduction process. 
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Figure 2: Upper and middle: Comparison of the total and inelastic neutron-nucleus cross sections calculated in the 
Glauber model with the available data. Bottom: The total pion- nucleus cross section as a function of the 

Glauber model calculations with the DL94 model for u^n are compared to the available data. 
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Figure 3: Comparison of the experimentally measured forward cross section of coherent p photoproduction on the 

n, lai, 1 11 1 


dll, |42i, |43f with the VDM-DL94 model and the Starlight parametrization. The red solid line 


proton 

shows the modified VMD (mVMD) parametrization (see text for details). 


This conclusion is confirmed by our observation that the latest 2006 HI data on the yp —)■ pp 
cross section 0 (we extrapolated the HI cross sections given at —t = 0.01 GeV 2 to — t = 0 
assuming the e Bt dependence with the value of the slope B reported by HI) disagrees with the 
normalization of the forward cross section calculated using the DL94 model by the factor of 0.84. 
This is seen in Fig. [3] where the forward yp —> pp cross section evaluated using Eqs. (E|) and (EJ) 
(the green dot-dashed curve labelled “VMD-DL94”) is compared to the whole bulk of the data. 
Also, for comparison, we show the parametrization of the forward yp —> pp cross section from 


the Starlight Monte-Carlo generator 48], which is widely used for predictions and modeling of 


vector meson photoproduction on nuclear targets. In order to agree with the 2006 HI data, the 
results of the VMD-DL94 and the Starlight parametrization should be decreased by the factor of 
approximately 0.7, which is much larger than what could be allowed by a variation of f p . From the 
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analysis presented above we can conclude the following: the assumption of the p meson dominance 
in the photon wave function has to be modified in order to agree to the whole set of data including 
the results of 2006 HI measurements. 

To this end, one can write the p meson photoproduction amplitude as the dispersion integral 
over the masses of the intermediate states generated in the 7 —> V transitions, which will involve 
the on-mass-shell fv, the pN cross section and the VN —» pN amplitude (here V denotes p-meson¬ 
like fluctuations of the photon with the invariant mass M, see our discussion in the Introduction). 
It is possible to demonstrate that inclusion of the contribution of the higher states can only weakly 
change f p , but it can noticeably reduce the cross section of the p meson production due non¬ 
diagonal transitions among different hadronic components of the photon and the p meson in the 
GVMD approach [fj, 10, 49j. On the other hand, within the VMD approach this can be modeled 
by defining the effective p-nucleon cross section a P N'. 


dpJV (Wiyp) 


fp A Q^. d a 'yp-+pp(t ~ 0) 

e V dt 


(9) 


We refer to this model as the modified vector meson dominance (mVMD) model; its prediction is 
shown by the solid red curve in Fig. [3] Note that a similar effect is also present in the CDM. 

The Gribov-Glauber model takes into account both elastic and inelastic diffraction; the latter 
leads to the additional—as compared to the Glauber model—inelastic nuclear shadowing contri¬ 
bution (the Gribov shadowing correction) [20J. The standard method to include this effect is 
given by the formalism of cross section fluctuations, which conveniently and successfully describes 
diffractive dissociation of protons, neutrons and pions on hydrogen and nuclei and inelastic nuclear 


shadowing in hadron-nucleus total cross sections 50 


Applying this formalism to the p meson-nucleus scattering, we obtain: 


mVMD-GGM 

<J ^A^pA 


- \ d 2 b 


daP(a) (1 - j 


( 10 ) 


which generalizes Eq. (]5j). 

The interpretation of Eq. (fTUj) is the following: the photon fluctuates into the p meson, which 
interacts with the target as a coherent superposition of eigenstates of the scattering operator, 
whose eigenvalues are the scattering cross sections a; the weight of a given fluctuation is given by 
the distribution P(er). Each state interacts with nucleons of the target nucleus according to the 
Gribov-Glauber model. The result is summed over all possible fluctuations, which corresponds to 
averaging with the distribution P(cr) at the amplitude level. 

Based on the similarity between the pion and p meson wave functions suggested by the additive 
quark model and our discussion above, it is natural to assume that P(cr) for the pN interaction 
should be similar to the pion P n (a), which we additionally multiply by the factor of l/(l + (cr/cro) 2 ) 
to take into account the enhanced contribution of small a in the pN interaction (we explained 
above that the contribution of small-cr fluctuations to the 7 IV —» pN amplitude is expected to be 
enhanced compared to the nN —> nN one): 


P(cr) = Cj 


3 -(o-/<T0-l) 2 /fF 


+ (^/ob) 2 


( 11 ) 


10 












The parameterization of Eq. (TIT]) satisfies the basic QCD constraint of P(a = 0) 7 ^ 0 and also 
P(a —> 00 ) —y 0. The free parameters C, ctq and hi are found from the following constraints: 


J daP(a) 


1 , 


J daP(cr)a = (a ), 

J daP(cr)a 2 = (tr) 2 (l + w„), 


( 12 ) 


where (a) = d p x in the mVMD model, see Eq. ((2]). 

The quantity c o a parametrizes the dispersion of P(cr) around its mean value (a), i.e., it charac¬ 
terizes the strength of cross section fluctuations. It can be determined using experimental informa¬ 
tion on the photon diffraction dissociation, in particular, the factorization of the photon and the 
pion diffraction dissociation cross sections scaled by the respective total cross sections. I 11 detail, 
the measurement [51] of inclusive diffraction dissociation of photons on hydrogen, yp —» Xp, in 
the range of 75 < E 1 < 148 GeV and M\/s <0.1 (M x denotes the produced diffractive mass) 
and the control measurement of inclusive diffraction dissociation of pions in the up —> Xp reaction 
at E n = 100 GeV showed that the respective M\ distributions scaled by the total cross sections 
are very similar in the photon and pion cases. For the cross sections integrated over M x , this 
observation means that: 


d(7-yp—}Xp(t 0 ) / dt dtj n p^xp(t 0 ) / dt uj 0 


a, 


ip 


<j. 


irp 


I 67 r 


&7tN 1 


(13) 


where in the last equation we expressed the cross section of pion diffraction dissociation in terms 
of uil characterizing the P 7 r (cr) distribution and the total pion-nucleon cross section cr nN . 

On the other hand, using the formalism of cross section fluctuations for the p-nucleon scattering 
and the mVMD model for the 7 — p transition, we obtain for the cross section of photon diffraction 
dissociation [compare to Eq. (15])]: 


d<r^p^xp(t 0 ) 

dt 


16vr \f 


2 r 


dcrP(a)cr 2 - (d pN f 


to, 


e 

16vr \ f p 


yVpN) 


(14) 


where the diffraction dissociation hnal state A" by construction does not contain p. The inelastic 
final state X is selected experimentally by analyzing the differential cross section as a function 
of the produced diffractive mass M x and corresponds to the values of M x beyond the p peak, 


M\ >1.5-2 GeV 2 


5l| . Substituting Eq. (TT4|) in Eq. (H3|) we obtain the desired constraint on u>, 

f 2 

J p ^ttN^'jp 


UJn = 




o 


(15) 


P N 


where the total photon-proton cross section a lp is taken from the fit to data [d|]. 

For the pion projectile, we use the constituent quark counting rule for the ratio of the nucleon- 
nucleon and the pion-nucleon total cross sections and obtain: 




2 u. 


N 


(16) 
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Here we effectively use validity of the limiting fragmentation which is well established experimen¬ 
tally. 

The pattern of cross section fluctuations for the nucleon projectile has the following dependence 
of the invariant collision energy y/s: the cross section fluctuations reach a broad maximum for 
24 < v'i < 200 GeV, are most likely small for y/s < 24 GeV and gradually decrease for y/s > 200 
GeV toward the Tevatron and LHC energies. Therefore, we use the following parametrization for 
the parameter c o^ describing the dispersion of the fluctuations: 


co 


N 

a 




PV i/24, 

P, 

P-0.15 ln(Vs/200) + 0.03(ln(y / s/200)) 2 , 


< 24 GeV, 

24 < y/s < 200 GeV, 
y/s > 200 GeV . 


(17) 


where the parameter f3 ~ 0.25 — 0.35 was determined from the analysis of pp and pp data 28 



20 40 60 80 


GeV 

Figure 4: The g 1 a^ p a cross section as a function of W 1V . The theoretical predictions using the mVMD model 
for the 7 p -A- pp cross section and the Gribov-Glauber model with cross section fluctuations for the 'yA -A pA 
amplitude are compared to the STAR (circle) and ALICE (triangle) data. The shaded area reflects the theoretical 
uncertainty associated with the parameter /3 characterizing the strength of cross section fluctuations (see text for 
details). 


It is known 22] from studies of corrections to the Glauber model for total proton-nucleus cross 
sections that suppression due to the inelastic shadowing is almost compensated by the effect of 
short-range correlations (SRC) in the wave function of the target nucleus. We included the effect 
of SRC by the following replacement 


T A (b) -> T A (b) + iAf- / dzp\(b, z ), 


(18) 
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where £ c = 0.74 fm is the correlation length. 

Onr predictions for the 7 A —> pA cross section as a function of W 1P are presented in Fig. [0 
The shaded area spanned by two red curves presents the results of the calculation using the 
mVMD model for the yp —> pp cross section and the Gribov-Glauber model with the effect of 
cross section fluctuations, see Eq. (ITU . The shaded area shows the uncertainty of our calculations 
due to the variation of the fluctuation strength oj a by changing /3 in the range 0.25 < (3 < 0.35. 
Our predictions are compared to the STAR (circle) and ALICE (triangle) data. One can clearly 
see from the figure that the inclusion of the inelastic nuclear shadowing enables us to explain 
the discrepancy between the UPC data on coherent p photoproduction on nuclei at large W 1P 
and the theoretical description of this process in the framework of the VMD-GM with the DL94 
parametrization of the pN cross section. 


4. Discussion 


The effect of the inelastic shadowing correction, which we demonstrate in these calculations, 
can be checked in the UPC measurements at the LHC. The inelastic nuclear shadowing changes 
the rapidity distribution of coherent p photoproduction in ion UPCs. Figure [5] presents the results 
of our calculation of do Pb p b ^. p p bPb /dy, see Eq. (jT() , as a function of the p meson rapidity y in Pb-Pb 
UPCs at the LHC at y/s^N = 2.76 TeV. The shaded area spanned by two red curves corresponds 
to the combination of the mVMD model and the Gribov-Glauber model for nuclear shadowing 
with cross section fluctuations (the shaded area shows the uncertainty of the calculations related 
to the variation of the fluctuation strength due to the change of /3 in the range 0.25 < /3 < 0.35); 
the blue dashed curve is the result of the calculation in mVMD-GM, i.e. without cross section 
fluctuations; the green dot-dashed curve is the result of the VMD-DL94 model combined with the 
Glauber model. The shape of the rapidity distribution predicted by the mVMD-GGM calculations 
is due to specifics of symmetric UPCs and the interplay between the energy dependence of the 
inelastic shadowing correction and the photon flux. 

The predicted shape of dop b p b ^ p p b p b /dy is different from the almost flat dop b p b ^ p p b p b /dy 
distribution obtained in the VDM-GM and Starlight approaches and is also in stark contrast 
with the calculations [53, [54] in the color dipole model approach predicting a bell-like shape for 
dop b p b ^ p p b p b /dy with the maximum at y = 0 and small values of dop b p b ^ p p b p b /dy at y ~ —4.5 
corresponding to W 1P ~ 5 — 10 GeV, i.e., to the energy range of the STAR measurements. From 
Fig. [4] it is seen that the experimental photoproduction cross section is almost constant in the 
energy range spanning the STAR and ALICE energies, Oryp b ^. p p b ~ 2 mb. In UPCs at y = 0, 
the contributions from both colliding nuclei serving as a target are equal, while at | y |= 4.5 the 
contribution of the low energy photon dominates. The photon fluxes are calculated in all studies 
similarly and with good accuracy, N^/ Pb (y = 0) = 108 and N^/ Pb (y = —4.5) = 250. Then one 
easily obtains that op b p b ^p b p bp (| y |= 4.5) ~ 500 mb > o Pb p b ^.p bPbp (y = 0) ~ 430 mb. These 
estimates confirm that the two-bumped shape of the rapidity distribution seems to be reasonable. 

The good agreement with the ALICE result allows us to predict the value of the cross section 
of coherent p photoproduction in Pb-Pb UPCs at y/sprN = 5.02 TeV in Run 2 at the LHC: 


da{y = 0 ) 
dy 


= 560 ± 25 mb . 


(19) 
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Figure 5: The rapidity distribution of coherent p photoproduction in Pb-Pb UPCs at tJsnn = 2.76 TeV. Theoretical 
predictions of the mVDM-GGM (red solid curves with the shaded area showing the uncertainty due to the variation 
of the fluctuation strength), the nrVMD-GM (blue dashed curve) and the VMD-GM (green dot-dashed curve) are 
compared to the ALICE data (see text for details). 


Examining the calculations of clastic photoproduction of p mesons on nuclei in the dipole model 


framework 53, [54], one notes that some of them describe the STAR and ALICE data while others 
do not — the results strongly depend on the models used for p-meson wave function and the dipole 
cross section. The dipole model framework was successfully used in the analyses of many processes 
studied at ffERA, such as, e.g., deep inelastic scattering (DIS) and vector meson electroproduction 
in a wide range of the photon virtualities Q 2 . ffowever, in processes dominated by soft physics such 
as, e.g., in photoproduction of light vector mesons, the application of the dipole approach is subject 
to signihcant theoretical uncertainties including the need to model the large-size contribution to 
the dipole cross section and the vector meson wave function. Considering the CDM predictions for 
p photoproduction in UPCs one finds that it is difficult to describe simultaneously the yp —>• pp 
and yA —> pA cross sections. Note also that the answer is sensitive to the assumed effective quark 
mass which enters in the photon wave function. Also, the use of a light quark mass (for example, 
~ 10 MeV in Qj) in several dipole models leads to a very large transverse size of the photon 
wave function and, consequently, to the f-dependence of the Compton clastic scattering which is 
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stronger than that observed in the data 55 . 

We also would like to briefly comment on the description of the STAR and ALICE data in the 
Starlight Monte-Carlo generator. From our point of view, the observed agreement with the data on 
coherent p meson photoproduction on nuclei is an effect of lucky coincidence and the weak energy 
dependence of the photoproduction cross section in the energy range covered by the STAR and 
ALICE measurements. The Starlight calculations are based on the parametrization of the forward 
IP PP cross section, VMD and the optical theorem. The 7 A —>• pA cross section is calculated 
in Starlight using the following expression: 


(f'yA—tpA 


d&'yA—tpA (t 0 ) 

dt 


F A (t) 1 2 dt = 


16 vr fl 


\<y 


totl 2 


F, i(f) 1 2 dt, 


( 20 ) 


where F A {t) is the nuclear form factor normalized by the condition Fa( 0) = 1. Note that the 
factorized form in Eq. (T20l) is an approximation. While the forward photoproduction cross section 
on a nuclear target in Eq. (T20l) follows from the VMD model, the optical theorem and the Glauber 
model (see Eq. (Q}), Eq. (l20]i does not take into account that the strong absorption of the p meson 
in the central region of a heavy nucleus results in narrowing of the momentum transfer distribution 
compared to that dictated by the undistorted nuclear form factor. I 11 particular, for heavy nuclei 
the first diffraction minimum is shifted by 10 — 15% compared to the position of the dip in the 
nuclear form factor. This shift is clearly revealed in the momentum transfer distributions at the 
rapidity y = 0 calculated in mVMD-GGM and Starlight approaches, which are shown in Fig. [61 
The effect of this shift appears to be qualitatively revealed by the STAR (Fig. 2 of 131) and ALICE 
(Fig. 3 of 0) results. 

A more serious shortcoming of the Starlight calculations is the use of the inelastic p-nucleus 
cross section instead of the total cross section in Eq. (l20li (compare the expressions in Eq. (( 8 ])), 
which violates the optical theorem. At high energies the cross section of elastic hadron scattering 
by a heavy nucleus is about 30% of the total cross section (see Fig. [ 2 ]). As a result, this decreases 
the estimate of the forward cross section by a factor of about two. One can see from Fig. [3] that the 
parametrization used in Starlight is close to that given by the DL94 model. Therefore, dividing 
the value of the VMD-GM cross section (green dot-dashed line) at y = 0 by the factor of two we 
reproduce the Starlight agreement with the ALICE data. 

Finally, we would like to emphasize that our calculations use as input the data on the forward 
p photoproduction cross section off the proton and the data on the photon and pion diffraction 
on hydrogen in a wide range of energies 5 GeV < W 1P < 50 GeV. As one can see from Fig. [31 the 
forward yp —» pp cross section is known in the the 5 — 20 GeV range with rather large experimental 
errors, while a p N extracted from these data is important for calculation of the cross section for the 
STAR kinematics and for predictions of the rapidity distribution at large | y | at the LHC energies. 
Some of these measurements would be doable with the recoil detector at COMPASS. Information 
at higher energies could be obtained from the studies of UPCs in pA at the LHC. It would be 
also very important to collect data on the photon and pion coherent diffraction off the proton and 
nuclear targets at high energies since only a handful of such data is available now. In the case of 
the yA process, one could get this information from UPCs at the LHC. 
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Figure 6: The momentum transfer distribution of coherent p photoproduction in Pb-Pb UPCs at y/S nn = 2.76 
TeV. The mVMD-GGM prediction (red solid curve) is compared to the Starlight result (blue dashed curve). 


Other possible directions of studies include coherent </> production, where we expect a significant 
amplification of the inelastic intermediate state effects due to the small value of cr(4>N) (such a 
measurement is certainly challenging for the main decay channel of (j) into two kaons, but the 15% 
pTT channel may work). 

The phenomenon of fluctuations of the interaction strength, which we discussed in the context 
of p meson exclusive production on nuclei, should also be manifested in a wide range of high energy 
7 A inelastic processes that could be studied in UPCs at the LHC. Effects of such fluctuations in 
inelastic pA collisions were considered in 56|] with experimental evidence reported in 


5. Conclusion 

With an increase of the collision energy, the composite structure of the photon becomes pro¬ 
gressively more pronounced, which leads to the following two features of the calculation of the cross 
section of p photoproduction on nuclei compared to the lower energies. First, the significant cross 
section of photon inelastic diffraction results in the sizable inelastic nuclear shadowing correction to 
the ykL —> pA cross section. Second, the QCD-motivated enhancement of the hadronic fluctuations 
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of the photon reduces the cross section of p photoproduction on the proton in agreement with the 
latest 2006 HI data. We took these features into account by combining the modified VMD model 
with the Gribov-Glauber model for nuclear shadowing, where the inelastic nuclear shadowing is 
included by means of cross section fluctuations described using a QCD-motivated parametrization. 
The resulting approach allows us to successfully describe the data on clastic p photoproduction on 
nuclei in heavy ion UPCs in the 7 < W 1P < 46 GeV energy range and to predict the value of the 
cross section of coherent p photoproduction in Pb-Pb UPCs at y 'snn = 5.02 TeV in Run 2 at the 
LHC: d<jpbPb^ P PbP b (y = 0 )/dy = 560 ± 25 mb. 
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